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4,  INTRODUCTION 

The  HBCU/MI  Partnership  Training  Award  project,  “  Photonic  Breast  Tomography  and 
Tumor  Aggressiveness  Assessment''’  is  designed  to  establish  a  breast  cancer  training  and  research 
program  at  the  City  College  of  New  York  (CCNY)  through  close  collaboration  with  the 

researchers  at  the  Memorial  Sloan  Kettering  Cancer  Center  (MSKCC).  The  focus  of  the  trainin  g 
component  of  the  project  is  to  familiarize  the  CCNY  researchers  who  happen  to  be  physical 
scientists  and  engineers  to  the  biological  aspects  of  cancer  research  through  attending  relevant 
courses,  and  cancer  research  practicum  thr  ough  laboratory  rotations.  The  objectives  of  the 
research  component  of  the  project  are  to  deve  lop  optical  imaging  and  spectroscopic  approaches 
to  (a)  distinguish  between  aggr  essive  and  slow  growing,  metast  atic  and  non-m  etastatic  tumors, 
(b)  non-invasively  detect  and  diagnose  breast  tumors  at  early  stages  of  growth. 

While  the  emphasis  of  the  program  is  on  the  training  com  ponent  during  the  first  two  years, 
significant  progress  has  been  m  ade  in  both  training  and  research  during  the  second  reporting 
period  (June  15,  2008  -  June  14,  2009)  covered  by  this  report. 

5.  BODY 

The  tasks  performed  and  the  progresses  made  during  the  current  reporting  period  are  as  follows: 

•  Accomplishment  of  pi  anned  training  objectives  through  laboratory  rotations  and 
attending  group  meetings  and  seminars  at  MSKCC;  and 

•  Pursuing  research  on  developm  ent  of  non-i  nvasive  optical  im  aging  a  nd  spectroscopic 
approaches  for  breast  tumor  detection. 

We  provide  a  brief  outlin  e  of  our  accom  plishments  in  th  ese  a  reas,  a  nd  refer  to  appended 
materials  for  detailed  description  where  applicable. 

5.1.  Accomplishment  of  Training  Objectives 

While  the  Yearl  training  involved  accumulation  of  basic  knowledge  of  the  biological  aspects  of 
cancer  research  through  course  work,  se  minars,  lectures,  and  workshops  [1],  the  e  mphasis  in  the 
second  year  was  on  hands-on  training  through  rota  tions  in  several  breast  cancer  research 
laboratories  at  MSKCC  {Specific  Aim  0,  Task  2  and  Task  3),  and  keeping  abreast  o  f  the  recent 
developments  in  breast  cancer  research  through  participation  in  group  m  eetings  and  sem  inars 
{Specific  Aim  0,  Task  5).  The  specialized  training  strategy  included:  training  on  animal  handling 
at  the  Research  Animal  Resource  Center  (RARC)  at  MSKCC,  and  rotations  through  laboratories 
engaged  in  sm  all  animal  imaging,  and  m  odern  bi  omedical  research  m  ethods.  Specifics  of  the 
training  received  are  described  below. 

5.1.1.  Animal  Handling  Training 

The  trainees  have  attended  the  following  small  animal  training  sessions  {Specific  Aim  0,  Task  2), 
held  by  RARC  at  MSKCC: 

•  RARC  Orientation  training  session; 

•  RARC  Basic  training  session  that  included  ha  ndling  and  re  straint;  inje  ctable 

anesthetics  and  anesthetic  m  onitoring;  parenteral  adm  inistration  (SC,  IP,  IV);  blood 
sampling;  and  euthanasia; 
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•  RARC  Rodent  Survival  Surgery  training  session  that  invol  ved  surgical  equipm  ent  and 
instruments,  disinfection  and  sterilization,  animal  preparation,  surgeon  preparation,  suture 
materials  and  wound  closure,  and  post-operative  care; 

•  RARC  Hazardous  Material  training  session  that  included  introduction  to  the  designated 
hazardous  m  aterials  use  areas,  entry  requireme  nts,  investigative  staff  responsibilities; 
animal/material  manipulations,  waste  disposal/transportation  of  materials,  and  euthanasia. 

•  RARC  Xenograft  training  session  that  focused  on  prevention  of  infectious  disease 
transmission  to  laboratory  re  searchers,  an  imal  care  s  taff  and  oth  er  indiv  iduals  who 
service  the  areas  where  these  experiments  are  performed; 

•  RARC  special  training  on  tail  vein  injection',  and 

•  RARC  special  training  on  necropsy. 

5.1.2.  Laboratory  Rotations 

The  trainees  rotated  through  th  e  following  four  laboratories  a  nd  received  training  on  im  aging 
techniques,  cell  culture,  bioluminescence  assay.  Western  Blotting,  and  other  key  techniques  used 
in  cancer  research. 

Rotation  I 

Laboratory  PI:  Dr.  Ronald  Blasberg 

Laboratory  background:  The  focus  of  the  current  research  in  this  lab  is  on  transg  ene  (repo  rter 
gene)  imaging  using  noninvasive  nuclear  and  optical  techniques. 

Training  Topic:  Cell  culture  and  Bioluminescence  Imaging 

•  Cell  culture  and  bioluminescence  training:  The  trainees  learned  th  e  basic  cell  culture 
techniques.  Each  of  them  had  the  responsibil  ity  to  m  aintain  two  genetically  m  odified 
breast  cancer  cell  lin  es:  MDA231-Fluc  -I-GFP  and  MDA468-Fluc-I-GFP.  Routine 
maintenance,  such  as  subculture  and  cell  counting,  was  performed. 

•  Bioluminescence  assay  and  imaging  training:  The  tra  inees  learned  the  lucif erase  assay 
method  for  evaluating  the  level  of  expres  sion  of  biolum  inescence  reporter  gene 
(luciferin)  and  learned  how  to  measure  the  linearity  of  bio  luminescence  intensity  cell 
numbers  using  Biospace  Photon  Im  ager.  One  of  the  tr  ainees  iden  tified  sa  turation 
problems  in  the  im  aging  system ,  which  m  ay  cause  m  isinterpretation  of  experim  ental 
results  for  large  num  her  of  cells  in  plates  or  animals.  He  suggest  ed  a  rem  edy  that  has 
been  implemented  to  correct  for  the  problem. 

Rotation  II 

Faboratory  PI:  Dr.  David  Sold 

Faboratory  Background:  This  labo  ratory  focuses  on  the  developm  ent  of  cancer  th  erapies  that 
target  pathways  responsible  for  cancer  initiation  and  progression. 

Training  Topics:  Western  Blotting  and  Flow  Cytometry 

The  techniques  that  the  trainees  learned  include: 
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•  Western  blot  technique  for  detecting  specific  p  roteins  rela  ted  to  RAS/RAF/MEK/ERk 
signal  pathway:  pERK,  ERk,  pMEK,  MEK,  Cy  clin  Dl,  p27,  RB,  PTEN,  Cleaved  PARP, 
and  Caspase-3. 

•  Drug  (PD901)  treatment  and  cell  growth  curve  m  easurement  for  0VCAR5,  T24,  and 
UMUC3  cell  lines. 

•  Flow  cytometry  technique  for  cell  cycle  analysis.  Nusse  nuclear  preparation  m  ethod  for 
DNA  content  was  taught. 

Rotation  HI 

Small- Animal  Imaging  Core 
Eacility  Manager:  Dr.  Pat  Zanzonico 

Eacility  Background  :  T  he  im  aging  core  facility  provide  s  MSKCC  investig  ators  with  unique 
capabilities  for  noninvasively  de  tecting,  localizing,  and  biological  ly  characterizing  prim  ary  and 
metastatic  cancer  cells  in  vivo  in  small-animal  models. 

Training  Topics:  Major  small  animal  imaging  techniques  and  systems 

The  trainees  learned  how  to  use  the  following  imaging  equipment  and  techniques: 

•  Micro-PET  PET  (positron  cm  ission  tom  ography  )scanners  from  CTI-Concorde 
Microsystems  ; 

•  Bioluminescence  imaging  system  from  Xenogen  (Capiper)  Corp.; 

•  Micro-CT  CT  (computed  tomography)  scanner  from  Imtek  Corp.; 

•  SPECT  /  CT  scanner  from  Gamma  Medica; 

•  Ultra  Sound  from  Visual  Sonics;  and 

•  Eluorescence  imaging  system  from  CRI. 

Rotation  IV 

Magnetic  Resonance  Imaging  (MRI)  Eaboratory  /  Small-animal  MRI  Core 
Eaboratory  PI:  Dr.  Jason  Koutcher 

Eaboratory/Eacility  Background  :  The  MRI/MRS  laboratory  /  f  acility  a  re  dir  ected  by  Dr. 
Koutcher  (Co-PI)  and  equipped  with  a  4.7  T  and  a  7  T  m  agnet  resonance  imaging  /spectroscopy 
systems. 

Training  Topics:  MRI / MRSI  techniques 

The  trainees  participated  in  the  following  procedures  and  experiments: 

•  Design  and  build  radiofrequency  coils  for  MRI/MRS  experiments; 

•  Acquisition  and  reconstruction  of  high  reso  lution  1-H  (standard)  im  ages  in  different 

orientations  (axial,  sagittal,  and  coronal)  in  sm  all  animals  (up  to  cat  size)  and  scans  on 
mice  whole-body,  brain,  lung,  kidney,  prostate,  bladder  cancers; 
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•  Dynamical  contrast  imaging  (DCE)  technique  for  analysis  of  perfusion; 

•  Magnetic  resonance  speetroseopy  of  lactate  for  metabolite  quantization  of  animals. 

5.1,3,  Other  Training  Activities 

The  CCNY  and  MSKCC  group  met  regularly  to  diseuss  the  progress  during  the  Year  2  training. 
The  trainees  attend  the  weekly  MSKCC  MRI  group  meeting  and  seleeted  seminars  {Specific  Aim 
0,  Tasks). 

5,2,  Development  of  Near-Infrared  Optical  Imaging  Modality  for  Breast  Cancer  Detection 

The  researeh  planned  to  be  pursued  in  the  project  has  two  components:  (a)  development  of 
non-invasive  near-infrared  optical  imaging  modalities  for  early  deteetion  of  breast  caneer,  and 
(b)  assessment  of  aggressiveness  of  tumor  growth  using  an  animal  model.  The  researeh  plan  was 
organized  such  that  work  on  developing  the  optical  imaging  modality  would  be  undertaken  early 
on,  while  the  work  on  animal  model  would  begin  in  the  third  year  of  the  project. 

Consequently,  the  work  on  development  of  non-invasive  near-infrared  optical  imaging 
modalities  for  early  detection  of  breast  cancer  {Specific  Aim  4)  that  started  during  Year  1  of  the 
projeet  was  further  pursued  throughout  the  eurrent  reporting  period.  This  researeh  builds  on  and 
extends  the  work  that  the  CCNY  group  has  been  pursuing. 

The  goal  of  the  research  is  to  develop  optieal  speetroseopy  and  imaging  approaches  that  use 
the  near-infrared  (800-1300  nm)  light  to  obtain  three-dimensional  (3-0)  tomographie  images  of 
human  breast  that  enable  detection,  localization,  and  possible  diagnosis  of  tumor(s)  in  the  breast. 
The  work  at  the  developmental  stage  would  be  carried  out  on  phantoms  that  have  optieal 
absorption,  emission,  and  scattering  properties  similar  to  those  of  breast  tissues,  and  on  realistic 
breast  models  assembled  using  ex  vivo  breast  tissues.  During  this  reporting  period  we  further 
extended  the  scope  and  explored  the  efficacy  of  the  optical  tomographic  imaging  using 
independent  component  analysis  (OPTICA)  [2-5]  that  we  developed  earlier,  and  explored  new 
approaehes  optical  imaging  using  principal  component  analysis  (OIPCA)  and  time  reversal 
optical  tomography  (TROT). 

5.2,1,  ‘Model  cancerous  breast’  studied  using  OPTICA 

We  extended  OPTICA  to  include  multi-wavelength  probing  and  explored  its  efficaey  on  a 
more  realistie  model  breast  {Specific  Aim  4,  Task  #15  and  Task  #16).  The  experimental 
arrangement  for  OPTICA,  which  we  assembled  {Specific  Aim  4,  Task  #14)  during  the  first 
reporting  period  [1],  uses  multi-source  illumination  of  sample  under  investigation,  and  multi- 
deteetor  transillumination  signal  acquisition.  Light  beams  of  wavelengths  750  nm,  800  nm,  and 
850  nm  from  a  Ti-sapphire  laser  were  used  for  probing.  The  ‘model  caneerous  breast’  was  a  100 
mm  diameter  and  42  mm  thick  cylinder  assembled  using  of  ex  vivo  adipose  tissues  {{Specific  Aim 
4,  Task  #15)  with  two  tumors  (invasive  duetal  earcinoma)  embedded  within.  A  16-bit  1024  x 
1024  pixels  CCD  eamera  was  used  to  reeord  the  signal.  The  sample  was  scanned  across  the  laser 
beam  in  a  16x26  x-y  array  of  grid  points  and  a  two-dimensional  transmission  image  was  recorded 
for  eaeh  position  for  each  wavelength  to  meet  the  multi-source  multi-deteetor  (eaeh  pixel  of  a 
CCD  eamera  being  a  deteetor  element)  imaging  arrangement  required  for  OPTICA.  The  data 
aequisition  time  is  less  than  8  minutes  for  a  16x26  sean  at  one  given  wavelength.  The  resulting 
data  was  analyzed  using  the  numerical  algorithm  of  OPTICA  [2-3].  The  details  of  the 
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experimental  arrangement,  data  aequisition,  proeessing  and  analysis  methods  are  similar  to  that 
presented  in  the  First  Annual  Report  [1]  and  a  published  paper  [5]. 

The  approaeh  provided  the  loeations  of  both  the  tumors  in  three  dimensions  with  high 
aeeuraey.  Multi-wavelength  measurements  enabled  better  diserimination  of  tumors  from  other 
eomponents.  A  baek-projeetion  algorithm  enabled  estimation  of  the  eross  seetion  of  the  tumors. 
The  sample  was  further  investigated  using  MRI  and  the  results  are  in  good  agreement  with  the 
OPTICA  estimates. 

A  paper  entitled,  ""Multi-wavelength  optical  tomography  using  independent  component 
analysis''’  based  on  these  results  has  been  accepted  for  presentation  at  the  NIH  Inter-Institute 
Workshop  on  Optical  Diagnostic  and  Biophotonic  Methods  from  Bench  to  Bedside  2009.  A  copy 
of  the  abstract  that  was  submitted  appears  as  Appendix  1.  A  manuscript  is  under  preparation  for 
publication  in  a  professional  journal. 

5,2,2.  Development  of  new  optical  imaging  approaches 

Fast  and  accurate  methods  are  needed  for  detection  and  localization  of  tumours  in  breast,  and 
for  detection  of  margins  during  surgical  removal  of  breast  tumours.  We  have  pursued  the 
adaptation  and  development  of  two  approaches  Optical  Imaging  using  Principal  Component 
Analysis  (OIPCA)  and  Time  Reversal  Optical  Tomography  (TROT)  and  tested  those  using 
experimental  data  on  absorptive  targets  embedded  in  a  breast  tissue  simulating  model  medium. 

Principal  component  analysis  (PCA)  is  a  multivariate  statistical  method  that  captures  the 
variance  in  a  data  set  in  terms  of  principal  components  and  enables  finding  patterns  in  data  of  high 
dimension  [6].  The  components  are  assumed  uncorrelated  as  opposed  to  the  similar  Independent 
Component  Analysis  (ICA)  approach  that  assume  the  components  to  be  independent.  In  our 
application  involving  detection  of  targets  in  a  turbid  medium,  light  propagation  from  different 
targets  to  the  source  plane  and  detector  plane  is  considered  to  be  uncorrelated  and  decomposed 
from  the  recorded  data.  The  locations  of  the  objects  are  retrieved  by  fitting  the  principal 
components  of  intensity  distributions  (PCID)  to  Green’s  functions. 

The  TROT  approach  [7,  8]  is  based  on  multi-static  data  and  vector  subspace  classification  to 
eigenvectors  of  a  round-trip  matrix.  In  optical  imaging  application,  a  response  matrix  represents 
the  transport  of  light  from  multiple  sources  through  a  turbid  medium  with  embedded  targets  to  an 
array  of  detectors.  The  response  matrix  is  constructed  from  the  experimental  data.  The  ‘round-trip 
(RT)  matrix’  is  constructed  by  multiplying  the  response  matrix  by  its  transpose  matrix  for 
continuous  wave  illumination  (by  adjoint  matrix  for  frequency-domain  case).  Mathematically,  the 
RT  matrix  is  equivalent  to  transfer  of  light  from  multiple  sources  through  a  turbid  medium  with 
embedded  targets  to  an  array  of  detectors,  and  back,  and  is  similar  to  the  time-reversal  matrix  used 
in  the  general  area  of  array  processing  for  acoustic  and  radar  time-reversal  imaging  [9]. 

Both  OIPCA  and  TROT  approaches  use  the  same  multi-source  probing  and  multi-detector  data 
acquisition  scheme  for  obtaining  multiple  angular  views  of  the  target,  as  OPTICA  does.  To  test  the 
efficacy  of  these  two  approaches,  we  carried  an  experiment  using  Intralipid- 10%  suspension  in 
water  as  a  model  medium,  and  two  absorptive  targets.  The  concentration  of  Intralipid- 10%  was 
adjusted  to  provide  a  transport  mean  free  path  1.43  mm  and  an  absorption  coefficient  pa  = 

0.003  mm’'  at  785  nm.  It  was  contained  in  a  250  mm  x  250  mm  x  50  mm  transparent  plastic 
container.  The  targets,  two  ~  10-mm  diameter  glass  spheres  filled  with  indocyanine  green  dye 
solution  in  water  (pa  =1.15  mm"'),  were  placed  at  20  mm  and  25  mm  away  from  the  entrance  face. 
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respectively.  The  sample  was  scanned  at  11x12  grid  points  with  a  step  size  of  5  mm.  The 
transmitted  signal  at  785  nm  from  the  exit  plane  was  recorded  by  a  CCD  camera.  From  each 
image,  a  53.4  mm  x53.4  mm  window  was  cropped  out  and  then  5x5  pixels  of  each  image  were 
binned  to  one.  Then  OIPCA  and  TROT  analysis  were  performed  separately.  A  schematic  diagram 
of  the  experimental  arrangement,  and  information  on  data  analysis  and  retrieval  of  target  positions 
appear  m  Appendix  2. 

Both  OIPCA  and  TROT  approaches  were  able  to  detect  and  obtain  the  positions  of  the  two 
targets.  The  retrieved  positions  of  the  targets  agreed  well  with  known  positions,  and  are  compared 
in  Table  I. 


Table  I;  Target  positions  determined  by  OIPCA  and  TROT 


Objects  ] 

Cnown  Positions 
{x,  y,  z)  (mm) 

OIPCA  retrieved  positions 
{x,  y,  z)  (mm) 

TROT  retrieved  positions 
{x,  y,  z)  (mm) 

1 

(46.5,  14.5,  20) 

(46.7,  14.7,  19.5) 

(46.7,  14.7,21) 

2 

(9.2,  44.5,  25) 

(7.6,  44,  24.8) 

(7.6,  43.2,  25) 

Our  future  work  will  include  exploring  the  efficacy  of  these  approaches  to  detect  and  locate 
scattering  targets,  as  well  as,  detecting  tumours  in  realistic  breast  models,  as  was  done  for  OPTICA 
(see  Section  5.2). 

6.  KEY  ACCOMPLISHMENTS 

•  The  key  training  accomplishm  ent  includes  successful  introduction  of  physical  scientists  and 
engineers  of  CCNY  research  team  to  cance  r  biology  research  m  ethods  through  laboratory 
rotations  at  selected  MSKCC  research  laboratories. 

•  Key  research  accom  plishments  include:  (a)  de  monstration  of  the  efficacy  of  Optical 
Tomography  using  Independent  Component  Analysis  (OPTICA)  approach  for  detection,  3-0 
localization,  and  cross  section  im  aging  of  tw  o  turn  ors  inside  a  realistic  breast  m  odel 
composed  of  excised  breast  tissues  ;  (b)  explorat  ion  of  the  efficacy  of  approaches  based  two 
new  algorithms  for  target  detection  and  localization. 

7.  REPORTABLE  OUTCOMES 
Conference  Presentations 

(1)  M.  Alrubaiee,  M.  Xu,  S.  K,  Gayen,  V.  Longo,  and  R.  R.  Alfano,  “Multi-wavelength 
optical  tomography  using  independent  component  analysis,’’’’  accepted  for  presen  tation 
(Presentation  Number:  NIH09-NIH01-53)  at  the  NIH  Inter-Institute  Workshop  on  Optical 
Diagnostic  and  Biophotonic  Methods  from  Bench  to  Bedside  2009,  October  1-2, 
Bethesda,  MD. 

(2)  Binlin  Wu,  M.  Alrubaiee,  W.  Cai,  M.  Xu  and  S.  K.  Gayen,  ^''Optical  imaging  of  objects  in 
turbid  media  using  Principal  Component  Analysis  and  Time  Reversal  Matrix  Methods,  ” 
accepted  for  presen  tation  at  the  A  nnual  Meeti  ng  of  the  Optical  Society  of  America, 
October  11-15,  2009,  San  Jose,  California. 
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8.  CONCLUSION 

The  work  carried  out  during  this  reporting  period;  (a)  con  tinues  the  training  of  CCNY  research 
team  in  biological  and  m  edical  aspects  of  labo  ratory  breast  cancer  res  earch;  and  (b)  shows  th  e 
potential  for  noninvasive  detection  and  three-dimensional  localization  of  a  tumor  within  a  breast 
with  s  ignificant  accuracy.  The  con  trast  is  based  on  th  e  differences  in  the  light  scattering  and 
absorption  characteristics  of  the  tumor  and  normal  breast  tissue. 

“No  What  Section” 

•  The  Nation  al  Cance  r  Institute  (NCI)  has  identified  the  developm  ent  of  im  aging 
methodologies  as  an  extraordinary  opportunity  for  advancement  in  cancer  research.  Since  the 
background  of  the  CCNY  tea  m  is  in  physical  sciences  and  engineering,  the  training  they 
received  would  provide  them  with  necessary  laboratory  b  ackground  in  the  biology  of  cancer 
research,  an  d  help  d  evelop  a  know  ledgeable  m  ultidisciplinary  research  force  in  th  e  fight 
against  breast  cancer. 

•  A  recent  study  involving  35,3 19  patients  underscores  the  influence  of  primary  tumor  location 
on  breast  cancer  prognosis,  and  m  akes  it  im  perative  that  breast  cancer  detection  modalities 
obtain  three  di  mensional  (3-Z))  location  of  the  turn  or  relative  to  the  axilla  [  10].  The  current 
work  is  an  important  development  in  obtaining  3-D  location  of  a  tumor  within  the  breast. 

•  The  study  of  model  cancerous  breast  assembled  using  ex  vivo  breast  tissues  is  important  and 
essential  for  the  next  step,  in  vivo  optical  breast  imaging  involving  volunteers. 
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Gayen,  '^Optical  imaging  of  objects  in  turbid  media  using  Principal  Component 
Analysis  and  Time  Reversal  Matrix  Methods,  ”  aceepted  for  presen  tation  at  th  e 
Annual  Meeting  of  the  Op  tieal  Soeiety  of  Am  eriea,  Oetober  11-15,  2009,  San 
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Appendix  1 

To  be  presented  at  NIH Inter-Institute  Workshop  on  Optical  Diagnostic  and  Biophotonic  Methods 

from  Bench  to  Bedside  2009 


Multi-wavelength  optical  tomography  using  independent  component  analysis 


M.  Alrubaiee',  M.  Xu^,  S.  K,  Gayen',  V.  Longo^,  and  R.  R.  Alfano' 


*  Institute  for  Ultrafast  Spectroscopy  and  Lasers,  Physics  Department 
The  City  College  of  New  York,  138th  Street  at  Convent  Avenue,  New  York,  New  York  1003 1 
^  Department  of  Physics,  Fairfield  University,  1073  North  Benson  Rd.,  CT  06824 
^  Zuckerman  Research  Center  of  Memorial  Sloan-Kettering  Cancer  Center  1275  York  Ave.,  NY,  NY  10021 


Multi-wavelength  optical  imaging  using  independent  component  analysis  (MW-OPTICA)  of  a  realistic  breast  model 
composed  of  ex  vivo  human  female  breast  tissue  is  presented.  The  sample  was  a  100  mm  diameter  and  42  mm  thick 
cylinder  formed  predominantly  of  adipose  tissues  with  two  tumors  (invasive  ductal  carcinoma)  embedded  within. 
The  experimental  arrangement  used  750  nm,  800  nm,  and  850  nm  light  from  a  Ti-sapphire  laser  to  illuminate  the 
sample,  and  a  16-bit  1024  x  1024  pixels  CCD  camera  to  record  the  signal.  The  sample  was  scanned  across  the  laser 
beam  in  a  16x26  x-y  array  of  grid  points  and  a  two-dimensional  transmission  image  was  recorded  for  each  position 
for  each  wavelength  to  meet  the  multi-source  multi-detector  imaging  arrangement  required  for  OPTICA.  The  data 
acquisition  time  is  less  than  8  minutes  for  a  16x26  scan  at  one  given  wavelength.  The  resulting  data  was  analyzed 
using  the  ICA  formalism. 

The  approach  provided  the  locations  of  both  the  tumors  in  three  dimensions  with  high  accuracy.  Multi-wavelength 
measurements  enabled  better  discrimination  of  tumors  from  other  components.  A  back-projection  algorithm 
enabled  estimation  of  the  cross  section  of  the  tumors.  The  sample  was  further  investigated  using  x-ray  computed 
tomography  (CT)  for  comparison  with  and  testing  the  efficacy  of  optical  approach.  MW-OPTICA  results  are  in  good 
agreement  with  CT  results. 

The  research  is  supported  in  part  by  US  Army  Medical  Research  and  Materiel  Command. 
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Appendix  2 


Optical  Imaging  of  Objects  in  Turbid  Media  using  Principal 
Component  Analysis  and  Time  Reversal  Matrix  Methods 

Binlin  Wu\  M.  Alrubaiee*,  W.  Cai‘,  M.  Xu^  and  S.  K.  Gayen' 

‘The  Institute  for  Ultrafast  Spectroscopy  and  Lasers,  Physics  Department 
The  City  College  and  the  Graduate  Centre  of  The  City  University  of  New  York,  New  York,  NY  10031,  USA 
^Physics  Department,  Fairfield  University,  Fairfield,  CT  06824,  USA 
Email:  hwu(a>,sci.ccnv.  cunv.  edit 

Abstract:  Principal  Component  Analy  sis  and  Time  Rev  ersal  Matrix  methods  were  used  to  develop 
approaches  for  i  maging  of  targ  ets  in  turbid  m  edia.  Th  e  efficacy  is  d  emonstrated  b  y  imaging  t  wo  targ  ets 
embedded  in  Intralipid- 10%  suspension  in  water. 

@2009  Optical  Society  of  America 

OCIS  codes:  (1 10.01 13)  Imaging  through  turbid  media;  (170.3880)  Medical  and  biological  imaging 

1.  Introduction 

Fast  and  accurate  methods  are  needed  for  detection  and  localization  of  targets  embedded  in  turbid  media.  In  this  paper, 
we  present  two  algorithms  based  on  Principal  Component  Analysis  (PCA)  [1]  and  Multiple  Signal  Classification 
(MUSIC)  to  the  eigenvectors  of  a  Time  Reversal  (TR)  matrix  [2-3]  for  obtaining  images  and  location  infonnation  of 
targets  in  turbid  media.  We  refer  to  these  approaches  as  Optical  Imaging  using  Principal  Component  Analysis  (OIPCA) 
and  Time  Reversal  Optical  Tomography  (TROT). 

2.  Formalism  and  Experimental  Arrangement 

A  multi-source  probing  and  multi-detector  signal  acquisition  scheme  to  obtain  multiple  angular  views  of  the  target(s) 
was  used  for  both  methods.  When  using  OIPCA,  the  leading  eigenvalues  provided  the  number  and  strengths  of  targets. 
Light  propagation  from  different  targets  to  the  source  plane  and  detector  plane  was  considered  to  be  uncorrelated  and 
decomposed  from  the  recorded  data.  The  location  of  the  objects  was  retrieved  by  fitting  the  principal  components  of 
intensity  distributions  (PCID)  to  Green’s  functions.  When  using  TROT,  a  Time  Reversal  (TR)  matrix  was  constmcted 
by  multiplying  the  response  matrix  by  its  adjoint  matrix  (transpose  for  continuous-wave  (CW)  illumination).  The  signal 
and  noise  subspaces  were  calculcated  and  separated  using  a  method  similar  to  L-curve  regularization.  The  eigenvectors 
with  leading  non-zero  eigenvalues  of  the  TR  matrix  correspond  to  the  embedded  objects,  which  are  orthogonal  to  the 
vectors  in  the  noise  subspace.  The  vector  subspace  method,  MUSIC,  along  with  Green’s  functions  calculated  from  an 
appropriate  forward  model,  was  then  used  to  determine  the  locations  of  the  embedded  objects. 


Fig.  1 .  A  schematic  diagram  of  the  experimental  aiTangement  used  for  imaging  objects  embedded  in  a  turbid  medium.  Inset  shows  the  2D 
array  in  the  input  plane  that  was  scanned  across  the  incident  laser  beam. 

The  sample  shown  schematically  in  Fig.  1,  was  a  250  mm  x  250  mm  x  50  mm  transparent  plastic  container  filled 
with  Intralipid- 10%  suspension  in  water,  and  two  absorptive  targets  were  embedded  in  it.  The  concentration  of 
Intralipid- 10%  was  adjusted  to  provide  a  transport  mean  free  path  l,~  1.43  mm  and  an  absorption  coefficient  =  0.003 
mm  ’  at  785  run.  The  targets,  two  ~  10-mm  diameter  glass  spheres  filled  with  indocyanine  green  dye  solution  in  water 
(fa  =  L15  mm’),  were  placed  at  20  mm  and  25  mm  away  from  the  entrance  face,  respectively.  The  sample  was 
scanned  at  11x12  grid  points  with  a  step  size  of  5  mm.  The  transmitted  signal  at  785  nm  from  the  exit  plane  was 
recorded  by  a  CCD  camera.  From  each  image,  a  53.4  mm  x53.4  mm  window  was  cropped  out  and  then  5x5  pixels  of 
each  image  were  binned  to  one.  Then  OIPCA  and  TROT  analysis  were  performed  separately. 
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Appendix  2 


Pseudo  Image 


3.  Results 

The  OIPCA  generated  PCIDs  of  target  1  on  the  deteetor  plane  and  souree  plane  are  shown  in  Figs.  2(a,e).  Fitting  of 
PCIDs  to  Green’s  funetions  are  shown  in  Fig.  2(b,d).  Similar  images  were  obtained  for  target  2. 
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Fig  2.  For  target  1,  OIPCA  generated  intensity  distributions  on  the  detector  and  source  planes  are  shown  in  (a)  and  (c),  respectively;  Green’s 
functions  fit  to  the  vertical  spatial  profiles  through  the  maxima  are  shown  in  (b)  and  (d),  respectively.  TROT  generated  cross-section  pseudo 
image  is  shown  in  (e)  and  pseudo  value  profiles  through  the  target  along  x,  y  and  z  directions  are  shown  in  (f). 


TROT  found  the  absorptive  objeets  and  generated  pseudo  images  using  pseudo  speetmm  from  MUSIC.  The 
pseudo  image  and  its  profde  for  target  1  are  shown  in  Fig.  2(e,f).  Similar  images  for  target  2  were  obtained. 
Loeations  of  the  objeets  found  by  OIPCA  and  TROT  are  listed  and  eompared  with  known  loeations  in  Table  I. 


Table  I:  Target  positions  determined  by  OIPCA  and  TROT 


Objects  1 

Lnown  Positions 
{x,  y,  z)  (mm) 

OIPCA  retrieved  positions 
(x,  y,  z)  (mm) 

TROT  retrieved  positions 
(x,  y,  z)  (mm) 

1 

(46.5,  14.5,  20) 

(46.7,  14.7,  19.5) 

(46.7,  14.7,21) 

2 

(9.2,  44.5,  25) 

(7.6,  44,  24.8) 

(7.6,  43.2,  25) 

4.  Summary 

Both  methods  deteeted  the  targets  and  provided  their  loeations  that  are  in  good  agreement  with  known  positions. 
Experiments  involving  tumors  in  model  breast  assembled  using  ex  vivo  breast  tissues  are  in  progress. 

The  research  is  supported  in  part  by  US  Army  Medical  Research  and  Materiel  Command. 
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